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FOREWORD 

This  report  was  prepared  by  the  Flight  Research  Laboratory 
under  an  internal  Project  on  Expenditure  Order  No.  R— 466— 1—4. 
It  was  decided  to  investigate  in  detail  the  problem  of  damping 
and  the  response  characteristics  of  systems  combining  inertial 
and  ground  velocity  indication.  The  broad  aspects  of  damping 
had  been  treated  previously  in  AF  Technical  Report  N  a  6045 
dated  October  1949  entitled  "Earth’s  Radius  Pendulum.”  This 
work  has  been  done  under  the  direction  of  Dr.  ).  E.  Clemens, 
Chief  of  Physics  Research  Branch,  Flight  Research  Laboratory. 


ABSTRACT 


This  report  deals  with  the  possibilities  of  combining  inertial  and 
ground  velocity  vector  information  in  order  to  obtain  a  system  with  im¬ 


proved  overall  response. 

The  classification  of  the  title  of  this  report  is  UNCLASSIFIED. 


PUBLICATION  REVIEW 


Manuscript  copy  of  this  report  has  been  reviewed  and  found  satis¬ 
factory  for  publication. 
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xChief,  Flight  Research  Laboratory 
Research  Division 
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INTRODUCTION 


Due  to  careful  scientific  research  and  skillful  design  certain  inertial  systems 
have  achieved  such  a  state  of  perfection  that  they  have  given  excellent  performance 
in  a  number  of  test  flights.  These  test  flights  have  proven  the  usefulness  of  those 
systems  for  long  range  navigation  of  airplanes  and  guided  missiles. 

Recently  it  has  been  claimed  that  another  system  based  on  the  doppler  radar  can 
give  the  ground  speed  and  the  drift  angle  under  flight  conditions  to  such  a  degree  of 
accuracy,  that  open  loop  integration  of  these  values  will  lead  to  a  position  indication 
regarded  as  good  enough  on  which  to  base  a  long  range  navigation  system. 

This  claim  led  some  of  the  sponsors  of  the  inertial  system  development  to  the 
consideration  it  might  be  worthwhile  to  use  the  doppler  radar  as  an  additional  source 
of  information  about  the  velocity  vector.  They  figured  that  the  platform  damping  could 
be  increased  as  the  velocity  vector  information  improved. 

It  is  the  aim  of  the  Chapters  I  and  II  of  this  report  to  deal  with  an  inertial  system 
with  additional  information  concerning  the  velocity  vector,  and  to  deduce  its  transfer 
functions.  Special  emphasis  is  placed  on  the  influence  of  the  errors  of  the  velocity 
vector  information.  However,  the  deductions  in  this  report  are  not  restricted  to  the  case 
that  the  information  concerning  the  ground  velocity  vector  is  delivered  by  a  doppler  radar; 
they  are  equally  valid  for  any  other  source  of  information  (e.g.  systems  based  on  optical 
or  infra-red  scanning  etc.). 

The  sponsors  of  the  doppler  radar  development  intend  to  use  a  simple  inertial  plat¬ 
form  with  accelerometer  in  order  to  improve  the  short  time  indication  which  is  otherwise 
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disturbed  by  noise*  This  possibility  is  dealt  with  in  Chapter  III. 

The  long  range  accuracy  of  the  combined  system  is  determined  in  the  first 
case  only  by  the  inertial  platform,  and  in  the  second  case  only  by  the  doppler 
radar.  It  is  evident  that  the  long  range  accuracy  of  the  combined  system  cannot 
be  better  than  the  long  range  accuracy  of  the  dominant  subsystem. 
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'  CHAPTER  1 

PROBLEM  OF  PLATFORM  DAMPING 

As  is  well  known  the  damping  of  an  inertial  platform  is  a  difficult  problem.  The 
platform  should  be  damped  against  the  true  vertical.  However,  without  additional  in¬ 
formation  the  apparent  vertical  is  the  only  term  available  on  the  platform  to  damp 
against.  It  is  evident  that  every  change  of  the  angle  between  the  true  and  the  appar¬ 
ent  vertical  produced  by  a  change  in  acceleration  (jerk)  will  disturb  the  platform 
through  the  damping.  Thus  one  is  faced  with  the  dilemma  of  compromising  between 
a  smaller  damping  term  than  desirable,  but  having  a  low  forced  error  disturbance,  or 
a  damping  term  of  the  proper  value  with  a  larger  forced  error  disturbance. 

The  ground  velocity  vector  is  the  additional  information  needed  to  compensate  for 

the  error  introduced  by  the  damping  term.  The  authors  proposed  to  use  as  compensation 

the  indicated  airspeed  for  the  lack  of  the  velocity  vector1.  Such  a  system  is  described 
2 

by  Hutzenlaub  .  It  is  evident  that  the  indicated  airspeed,  representing  only  the  magni¬ 
tude  of  the  ground  velocity  vector  falsified  by  the  wind,  is  an  imperfect  substitute.  The 
recent  development  of  devices  to  measure  the  true  ground  velocity  vector  inspired  the 
authors  to  this  investigation. 
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CHAPTER  II 


SPIRE  SYSTEM  WITH  GROUND  VELOCITY  VECTOR  METER 
The  MIT— Spire  system  (see  Bibliography  No.  2)  shall  be  taken  as  an  example  of 
an  inertial  system.  This  system  is  aligned  with  respect  to  the  preselected  great  circle 
course  and  the  surface  of  the  earth,  and  it  is  automatically  controlled  so  that  it  retains 
this  alignment  all  the  time  within  very  close  limits.  The  sensing  elements  are  two 
accelerometers;  the  one  measuring  the  component  tangential  to  the  great  circle  course 
furnishes  the  input  for  the  range  computer,  the  other  one  measuring  the  component 
normal  to  the  great  circle  course  furnishes  the  input  for  the  tracking  computer.  Thus 
it  is  evident  that  the  tangential  component  of  the  velocity  vector  must  be  used  to 
compensate  for  the  disturbance  introduced  by  the  damping  of  the  range  computer,  and 
the  normal  component  of  the  velocity  vector  must  be  used  to  compensate  for  the  dis¬ 
turbance  introduced  by  the  damping  of  the  tracking  computer. 

Because  the  angle  A between  the  programmed  great  circle  course  and 
the  path  (see  Fig.  2b)  is  always  small,  the  tangential  component  is  the  true  ground 
speed  ^  A/C)1  anc*  t^ie  normal  component  is 

^T(E  -  A/C)  =  V(E  -  A/C):  A  (p  .  path)  (1) 

However,  the  angle  A  *  _  pa^  is  not  measured  directly,  since  the  ground  velocity 

vector  meters  give  the  ground  speed  ^  A/C)  and  the  an£le  (x  .  path)  ° *  the 
path  against  the  longitudinal  axis  X  of  the  airplane. 

To  obtain  the  desired  angle  path)  anS^c  A  *  (p  *  x)  ° *  t^ie  longitudinal 

axis  against  the  programmed  great  circle  course  which  is  always  indicated  without  error 
on  the  Spire  gimbal  system  must  be  added  to  ^  path)  (see  Fig.  2b). 
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A 


(P  •  X) 


(x  -  path) 


(2) 


=  A* 


A* 


(P  •  path ) 

The  velocity  components  (M)  A/C)  aDc^  W)  ^T(E  -  A/C )  are  furnished  by  the 

ground  velocity  vector  meter  obviously  have  some  error.  It  is  assumed  that  the  indicated 
values  can  be  represented  as  the  sum  of  a  correct  term  and  an  error  term  (EM)* 


Thus  one  obtains  the  equations 

(M)  V (E  _  A/C)  =  \(E  m  A/C)  +  (EM)  V (E  .  A/rC)  (3) 

and 

(M)  VTfB  _  A/C)  =  VTfB  .  A/C)  +  (EM)  VTfE  „  A/C)  (4) 

=  V  *A*  + 

y  (E  -  /1/C)  M  (P  -  p«*i>)  + 

+  ^  (P  -  pfl/ib)  '  ^(E  -  A/C)  + 

+  ^ (E  -  A/C)  '  (^)  \ x  -  path) 

The  term  with  the  product  of  the  two  errors  is  neglected  in  (4),  because  it  is  a  small 
second  order  term. 

1.  Transfer  Function  of  the  Range  Indicating  System  with  Ground  Velocity  Vector  Meter 
Fig.  la  shows  the  block  diagram  of  the  system  with  a  damping  loop  and  a  ground 
velocity  vector  meter;  which  corresponds  with  the  system  Fig.  V— A  on  page  143  of 
Hutzenlaub’s  report2. 

The  external  information  about  the  velocity  shall  be  disregarded  for  a  moment.  It 
is  well  known  that  such  a  platform,  if  it  has  no  damping  loop,  is  on  the  borderline  of 
instability,  but  it  is  insensitive  against  external  disturbances  (i.e.  against  changes  in 
acceleration).  On  the  other  hand,  such  a  platform  with  a  damping  loop  is  sensitive 
against  changes  in  acceleration. 
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The  disturbing  damping  loop  consists  in  the  feed  back  of  the  velocity.  If 
therefore,  as  shown  in  Fig.  la,  mn  additional  external  source  of  information  about 
the  velocity  is  used  to  compensate  for  the  disturbing  internal  velocity  feed  back, 
the  platform  behaves  like  an  undamped  one  with  respect  to  external  disturbances, 
but  nevertheless  any  oscillations  of  the  platform  (however  they  may  have  been 
originated)  die  out. 

The  transfer  function  can  be  derived  from  Fig.  la.  One  obtains 


_ 


^  -  s3  +  SoS2  +  s  +  S5 


(P  -  i) 


SiS3  +  S2s2  +  s  +  S5 


A  + 
m(p  -  o + 


s  -  Sa 


2  W  J 


(EM)  V 


SlS3  +  S2s2  +  s  +  S: 


m  -  a/ci 


(s) 


(5) 


The  terms  caused  by  the  correct  part  _  A/C)  ~  ^(p  -  t)  oi  t^le  grouod  velocity 

information  are  contained  in  the  first  expression  of  the  right  side  of  equation  (5).  The 
second  expression  comprises  only  the  influence  of  the  ground  velocity  meter  error. 

If  the  platform  is  tuned  to  84.4  minutes,  i.e.  if  one  has  made 

the  transfer  function  of  the  first  term  is  equal  to  one.  The  error  of  the  indication  may 
be  defined  in  the  usual  manner. 


(E)A 


=  A  “A 

(P  -  i)  (P  -  i)  (P  •  t) 


Thus  one  obtains  by  combining  (7)  with  (5)  for  the  tuned  platform  (Si  =  1/Wg) 


(6) 


(7) 
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:'»•  *  -p~ > 5 

©V-« - - 


(EM)  V 


1  3  #»  2  .  » 
J —  s  +  S2s  +  s  +  S2 


w 


ffi  -  A/C) 


(s) 


Finally  we  introduce  the  dimensionless  operator 


D  = 

W, 


(9) 


and  divide  by  .  A/c/^E’  °btaining 


Rfi(E)  V  -  «  1  {S’~  -yr]  WF  • 1 D  <EM)  V(E  -  A/c/D) 


(E  -  A/C) 


WP  D3+ S2WpD2+ D  + -Sa 

wB 


no; 


(H  -  A/O 


Comparing  the  transfer  function  (10)  with  the  transfer  function  of  the  system  without 
velocity  meter  information  as  it  is  written  on  page  154  of  Hutzenlaub’s  report,  one 
recognizes  that  the  error  is  reduced  to  the  value  of  the  product  of  this  error  and  the 
relative  velocity  error.  That  is  indeed  a  worthwhile  improvement. 

It  may  be  noted  that  the  transfer  function  equation  (10)  is  positive,  while  the 
transfer  function  for  the  system  without  velocity  information  is  negative.  This  is 
evident  when  one  bears  in  mind  that  the  velocity  meter  information  is  added, while 
the  feed  back  of  the  speed  $2S  A^p  ^  ^  is  subtracted  (see  Fig.  la). 

2.  Transfer  Function  of  the  Tracking  System  with  Ground  Velocity  Vector  Meter 

Fig.  2a  shows  the  block  diagram  of  the  system  and  Fig.  2b  the  nomenclature 
and  kinematics.  Investigating  the  tracking  system,  it  is  evident  that  the  velocities, 
angles  etc.  are  rectangular  to  the  plane  of  the  programmed  great  circle  course.  In 
the  formulas  this  fact  is  expressed  by  the  subscript  T.  (Concerning  w  A/C) 

see  the  formulas  (1)  to  (4)  ). 

The  transfer  function  of  the  tracking  computer  with  Schuler  tuning  (84.4  minutes 
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period)  for  a  cross  wind  disturbance  can  be  derived  from  Fig.  2a. 


R/1Aj:(p..z  *)-  =—  .  Fi(D)  -  — L  F2(D)  .  F3(D)  .  (gjj)  Vyre  ■  a/c)  M 

^T(E  •  air) 


V T(E  ■  air) 


with 


Fi 


*e 

V 

D 

(C3  +  1)D2+  - 

Cl  D  +  Ca 

WE  W2 

c3d2  +  Cl 

D  +  C‘ 

We 

We 

(C3+  1)D2  + 

Cl  D  +  Cz 

v/E  w£ 

(S2-  s> 

K 

>V> 

3  .  r  hi  f>2  ,  n  ,  ^3 


D  +  S2W£D  +  D  + 


W, 


The  first  term  of  the  right  side  of  (11)  is  independent  of  the  velocity  vector  meter 
error;  Fi  represents  the  transfer  function  of  the  control  computer.  The  character¬ 
istic  equation  (i.e.  the  denominator)  describes  a  second  order  system  with  the  mass 
term  (C3  +  1),  the  damping  term  Ci/W^  and  the  spring  term  C2/Wg  ;  but  one  must 
have  in  mind  that  the  equation  is  written  in  units  of  the  dimensionless  operator 
D  =  —  (which  is  the  ratio  between  the.  operator  s  of  the  system  and  the  Schuler 

frequency  =  4.47  hrs.’  ).  In  an  earlier  investigation  the  authors  found  that  a 
tracking  computer  without  ground  velocity  vector  information  gives  a  response  as 


(U) 
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good  as  possible  with  C3  =  10  to  20,  Ci  -  20.7  hrs  ,  C2  ~  24.8  hrs  2.  The  dimension¬ 
less  natural  frequency  was  0.22  and  the  damping  ratio  0.3* ** 

The  situation  is  much  more  favorable  for  the  system  with  ground  velocity  vector 
meter  because  of  the  achieved  decoupling.  C3  can  be  made  zero,  and  from  the  view¬ 
point  of  stability  there  is  no  objection  to  change  C*  or  C2  *n  order  to  obtain  a  faster 
response  and  a  better  damping  ratio.  With  Q±  =  20.7  hrs."1  and  C2  =  218  hrs."  one 
obtains  the  dimensionless  frequency  2.12,  and  the  damping  ratio  0.7. 

The  second  term  of  equation  (11)  represents  the  disturbance  of  the  system  due  to 
the  error  of  the  ground  velocity  vector  meter.  It  is  represented  by  the  product  of  two 
transfer  functions  F2(D)  and  F  3(D).  F2(D)  contains  the  constants  of  the  control  com¬ 
puter  and  F 3(D)  those  of  the  indication  computer.  F2(D)  approaches  1  if  C3  »1*  Then 
the  characteristic  equation  of  the  second  term  is  of  the  third  order.  The  condition 
for  stability  due  to  Routh’s  criterion  is 

S2Wg>S3  (12) 

This  can  be  fulfilled  eas;ly,  because  in  the  system  with  ground  velocity  vector  meter 
there  are  no  restrictions  to  the  damping  constant  S2. 

If  C  3  is  chosen  small  (e.g.  zero)  the  transfer  function  F2(D)  is  unequal  to  1.  Then 
the  characteristic  equation  can  be  represented  by  the  product  of  a  quadratic  and  cubic 
equation.  Because  of  the  factorization  the  stability  check  for  the  fifth  order  system  is 

*  The  authors  do  not  see  any  reason  why  almost  9  fold  Increase  of  the  gain^  should  lead  Into  difficulties, 
because  the  noise  is  very  low  in  this  loop  due  to  the  manifold  Integration.  Anyway,  a  large  increase  should 
at  least  be  possible. 

**  The  loop  with  S3  was  added  to  obtain  the  correct  response  on  the  low  frequency  side. 
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simple;  it  is  sufficient  to  show  that  both  factors  represent  stable  systems.  Besides 
the  requirement  that  all  coefficients  must  be  positive,  the  only  additional  requirement 
is  the  fulfillment  of  Routh’s  criterion  for  the  cubic  equation,  represented  in  (18). 

This  shows  that  the  tracking  computer  with  ground  velocity  vector  meter  can  be  made 
stable  very  easily. 

3.  Comparison  Between  the  Tracking  System  Without  and  With  Ground  Velocity  Vector  Meter 

At  this  point  it  might  be  interesting  to  compare  the  system  without  and  the  one 
with  ground  velocity  vector  meter  in  order  to  appreciate  fully  the  enormous  gain  in 

stability  of  the  tracking  computer  with  ground  velocity  vector  meter. 

2 

The  earlier  investigations  conducted  at  M.  I.  T.  had  revealed  that  a  tracking  system 
without  ground  velocity  vector  information,  and  with  indication  computer  tuned  to  84.4 
minutes  becomes  dynamically  unstable  the  more  the  gain  S2  of  the  damping  loop  is  in- 
creased.  The  authors  showed  that  stability  could  be  achieved  by  adding  a  damping  loop 
C3  to  the  control  computer,  but  even  in  this  case  the  response  is  not  very  good. 

In  Table  la,  page  9?  are  contrasted  the  coefficients  of  the  characteristic  equation  of 
the  system  without  and  the  one  with  ground  velocity  vector.  Those  coefficients  which  do 
not  coincide  in  both  cases  are  framed.  The  stability  of  the  system  can  be  deduced  from 
the  coefficients  in  Table  la,  page  9>  by  utilizing  Routh’s  stability  criterion.  However,  in 
the  case  of  the  fifth  order  systems  this  task  is  somewhat  cumbersome.  Therefore  the 
expressions  of  Routh’s  stability  criterion  are  written  down  for  the  simplified  case  without 
S3  and  C3  in  Table  lb,  page  9*  This  advantage  of  the  system  with  doppler  radar  is 
illustrated  by  Fig.  3. 
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CHARACTERISTIC  EQUATIONS  FOR  TRACKING  SYSTEM 


CHAPTER  HI 


DOPPLER  RADAR  WITH  INERTIAL  INFORMATION 
TO  IMPROVE  THE  SHORT  TIME  INDICATION 

As  will  be  pointed  out  in  Chapter  IV,  the  short  time  indication  of  the  doppler  radar 
is  disturbed  by  noise.  One  can  filter  out  the  noise,  but  then  the  response  of  the  system 
becomes  sluggish,  as  can  be  seen  from  Figure  4b,  curve  2  which  shows  the  response  to 
a  step  input  in  acceleration.  But  if  the  area  represented  by  the  difference  between  the 
curves  1  and  2  of  Figure  4b  could  be  added,  the  response  of  this  improved  system  would 
be  correct.  Indeed,  the  necessary  quantity  (see  curve  3  of  Figure  4b)  can  be  formed  from 
the  indication  of  an  accelerometer  mounted  on  a  horizontal  platform. 

The  block  diagram  of  this  system  is  shown  in  Figure  5«  The  first  order  filter  smooths 
the  doppler  velocity  indication;  the  lag  is  compensated  for  by  the  accelerometer  feed  in. 
The  transfer  function  can  be  easily  deduced  from  Figure  5: 

V.  =  JL  {CV  +  C(EDo)Y  -CV.  +  a  +  (E)  a  }  (13) 

*  s  * 

Having  in  mind  that  a  =  sV  one  obtains: 

V.  =  V  +  — %r-  (EDo) V  +  — L  (E)  a  (14) 

1  s  +  C  s  +  C 

This  equation  shows  that  the  indication  is  correct  at  every  instant,  except  for  the  errors 
introduced  by  the  doppler  radar  and  the  indication  of  the  acceleration.  It  can  be  shown 
that  the  latter  error  is  essentially  generated  by  the  misalignment  of  the  platform  on  which 
the  accelerometers  are  fixed.  If  the  angle  A(t-i)  between  the  true  vertical  and  the  vert¬ 
ical  indicated  by  the  platform  is  small,  one  obtains: 

(E)o  =  gA  (fmi)  (15) 
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In  the  previous  paragraph  the  platform  with  the  accelerometers  was  assumed  to  be  held 
horizontally,  but  no  consideration  was  given  to  the  fact  how  this  could  be  achieved,  In  this 
paragraph  shall  be  shown  how  the  combined  system  consisting  of  the  doppler  radar  and  the 
platform  can  be  linked  together. 

The  sponsors  of  the  doppler  radar  suggested,  if  one  could  use  a  simple  first  order  plat¬ 
form  (similar  e.g.  to  the  normal  airplane  horizon)  as  a  base  for  the  accelerometers,  since  a 
complete  second  order  platform  tuned  to  84.4  min.  means  a  very  high  expenditure  for  the 
purpose  in  question.  As  is  well  known  the  first  order  platform  develops  two  kinds  of  errors: 
it  tries  to  adjust  itself  into  the  apparent  vertical  during  periods  of  accelerated  flight,  and 
it  develops  misalignment  as  a  consequence  of  the  velocity  of  the  vehicle*.  The  latter  error 
orignates  from  the  fact  that  even  during  a  uniform  motion  along  the  surface  of  the  earth  a 
misalignment  of  the  platform  must  develop  in  order  to  generate  the  necessary  error  signal 
which  makes  the  platform  turn  with  the  rate  of  the  true  vertical.  Since  the  doppler  radar 
gives  information  about  the  velocity  vector  one  has  in  the  combined  system  the  possibility 
to  compensate  for  the  velocity  error  of  the  first  order  platform, 

2.  Transfer  Functions  of  the  Doppler  Radar  with  a  First  Order  Platform 

Figure  6a  shows  the  schematic  and  Figure  6b  the  block  diagram  of  the  combined  system. 
The  lower  part  of  Figure  6b  representing  the  doppler  radar  with  low  pass  filter  is  identical  with 
Figure  5*  The  upper  part  of  Figure  6b  is  the  complete  representation  of  the  first  order  plat¬ 
form  with  accelerometer  and  additional  velocity  feed  in. 

Two  transfer  functions  of  the  system  are  of  special  interest,  the  transfer  function  which 
relates  the  indicated  velocity  to  the  true  one,  and  the  transfer  function  which  relates  the  mis¬ 
alignment  of  the  platform  to  the  angular  acceleration. 

For  the  velocity  transfer  function  one  reads  from  Figure  6b  the  following  four  equations. 
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(16) 


is  +  CJ  V,  =  CXV  +  C ! (EDo) V  +  K2Rb k(i .  a)  +  K2R£  (E)  A 


^(i  -  A(o  -  a)  ~  A(o  -  *; 


A  -  1  .  s%  WH 

f° ' a>  R£s  W2 


a  7; 

as> 


A 


fo 


i) 


~^(i  -  +  -Me)  A  +  — C^_V. 

s  s  R£s 


(EDo)  Vi  s  the  errot  of  the  doppler  radar  and  (  E )  A  is  an  angular  error  in  the  alignment  of  the 
platform  which  might  be  superimposed  to  the  misalignment  of  the  platform  due  to  the  dynamics. 
This  additional  error  (E)  A  will  be  caused  e.g.  by  Coulomb  friction  in  the  gimbal  system  etc., 
thus  in  most  cases  it  will  be  a  region  of  uncertainty. 

This  uncertainty  in  the  platform  alignment  causes  an  uncertainty  in  the  indication  of  the 
acceleration.  Both  are  related  by  the  equation: 


(E)A  =  ■(  (20) 

9 

Combining  the  equations  (16)  to  (19)  leads  to: 

s2  +  CiS  +  (CjKi  +  K2)  C^S  +  KJ  (EDo)  V  +  K2  R£s (E)A 

V,  =  J- - V+  -  (21 

s2  +  (C1  +  )s  +  (CiKi  +  C2K2)  s2  +  (C±  +  Kx)  s  +  (CJ^  +  C2K2) 


The  transfer  function  which  relates  V .  to  V  should  become  unity.  First  of  all,  this  condition 
must  be  fulfilled  for  the  steady  state  condition  is 


C2  *  1  (22) 

The  additional  requirement  to  make  the  transfer  function  equal  to  one  on  the  very  high  end  of 
the  frequency  band  (s  00  )  is  evidently 
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(23) 


But  that  means,  as  shall  be  shown  later,  in  Chapter  E3.3,  the  tuning  of  the  whole  system 
to  the  Schuler  frequency. 

One  can  make  the  transfer  function  equal,  to  unity  for  all  values  of  S  by  fulfilling  the 
additional  condition 


Kx  =  0 


(24) 


Also  the  meaning  of  this  will  become  clear  later. 

Now  the  transfer  function  which  relates  the  misalignment  A^  w  ^  of  the  platform  to  the 
angular  acceleration  sV/Rj?  shall  be  deduced. 

A  =  A  _  A 

*(t  .  i)  R(o  -  i)  -  *(o  -  t) 

Using  the  relation 


A 


(O  .  t) 


(26) 


in  addition  to  the  equations  (17),  (18),  (19),  (21)  and  to  the  steady  state  condition  (22),  one  can 
deduce 


1  .  KjS  +  tCxKr+K^yg  sV  +  C1(EDo)V/Rg  +  {KxS+^Kx+K^KEjA 

Wg  sMCx+KJs  +  fCJVI^)  Rg  s^fCi  +  KJs  +  ^Kx+KJ 


(27) 


For  the  tuned  system  (i.e.  K2  =  Wg  )  one  obtains 

A*  1  Kxs  +  CxKx  ,  sV  +  C1(EDo)V/Rg+  (KxS  +  tCxKx  +  W*)  }(E)A  (28) 

A(t-i)~~  #  ’  * 

Wg  s^tCx+K^s+^xKx  +  Wg)  Rg  sMCx+KJs+fCxKx  +  W^) 

And  if  finally  is  made  K*  =  0,  the  first  right  hand  term  vanishes  completely,,  and' the  misalignment 
is  caused  only  by  the  errors. 
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(29) 


A 


** 

(t  -  i) 


C1(EDo)V/Rh  +  Wg  (E)  a 

s2  +  ClS  +  Wg 


In  conclusion  the  free  oscillation  of  the  system  shall  be  studied.  It  is  described  by  the 
denominator  of  equation  (27),  or  of  equation  (21)  and  considering  the  steady  state  condition 
(22).  The  general  characteristic  equation  of  a  second  order  system  is 


sz  +  2RW  s  +  w; 


(30) 


Wo  means  the  frequency  of  the  undamped  system,  and  R  is  the  damping  ratio.  As  well  known 
R  <1  characterizes  a  periodic  motion,  and  R  J>  1  a  periodic  one.  Comparing  the  coefficients 


of  the  denominator  of  (27)  or  of  (21)  with  (22)  with  those  of  (30)  one  finds 

n  1  C±  +  Kx 


2  \CjKi  +  Ks 


(31) 


W, 


C1K1  +  Ks 


(32) 


Therefore,  the  aperiodic  case  will  exist,  if 

4K2<(C1-K1)2  (33) 

or  for  the  tuned  system  (K2  =  Wg,  see  equ.  (23)X  if 

2Wb£  |C1-K1|  (34) 

In  order  to  obtain  a  good  transient  response  K 1  should  be  made  small  compared  with  C].. 

Since  on  the  other  hand  the  constant  Ci  of  the  doppler  radar  filter  is  very  large  compared  with 
wE,  the  combined  system  will  be  always  very  strongly  overdamped,  if  it  is  designed  to  have 
a  good  transient  response. 

Fig.  7  illistrates  the  transient  response  of  the  system  for  a  step  input  in  velocity.  As 
can  be  readily  seen  from  these  curves,  the  transient  response  can  be  improved  if  Kx  is  chosen 
small,  i.e.  if  the  time  constant  of  the  inertial  platform  is  long  as  compared  to  the  time  constant 

M 
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of  the  doppler  filter.  This  requirement  calls  for  high  quality  inertial  components. 

The  curves  in  Fig.  7  refer  to  the  principle  part  of  the  transfer  function  (21), 

describing  the  dynamic  response  of  the  system.  The  effects  of  the  errors  in  ground 

velocity  and  acceleration  measurements  (edo)  and  (e)d,  which  are  statistical  in 

nature,  have  to  be  superimposed  upon  these  response  curves.  The  RMS  —  error  in 

V  j  d  due  to  doppler  error  (sdo)V  with  the  filtering  time  used  in  the  examples  of 

Fig.  7(T  =  20  Sec)  is  less  than  0.1  %  (this  statement  is  based  on  the  results  of 

F.  B.  Berger,  reference  4).  The  RMS  —  error  in  V  .  ^  due  to  acceleration  (e)o 

will  depend  on  the  quality  of  the  inertial  components  used. 

3.  The  Combined  System  from  Another  Aspect 

Now  the  reason  for  the  behaviour  of  the  combined  system,  as  it  was  represented 

by  equations  (21),  and  (27)  to  (29),  shall  be  investigated.  For  this  purpose  the  system 

♦ 

is  shown  in  Figure  8  in  two  different  representations.  Figure  8a  is  an  identical  replot 
of  Figure  6b,  while  Figure  8b  shows  the  same  block  diagram  in  another  arrangement. 

As  one  can  see  all  the  loops  are  completely  alike  in  both  figures.  But  Figure  8b  re¬ 
veals  at  once  that  this  combination  of  the  two  first  order  systems  (which  was  proposed 
as  a  simple  arrangement  in  the  beginning  of  this  paragraph)  is  essentially  identical 
with  a  damped  second  order  platform  with  doppler  velocity  information.*  It  is  evident 
that  this  second  order  platform  can  give  the  correct  response  during  a  transient  state 
only,  if  tuned  to  the  Schuler  frequency.  One  recognizes  further  that  the  system  has 
two  damping  loops,  a  feed  back  loop  with  the  gain  C±  and  a  forward  loop  with  the  gain 
Kj.  The  disturbance  of  the  system  by  the  feed  back  loop  is  compensated  by  the  doppler 

*  This  results  Iron  the  fact  that  two  crossfeeds  have  been  introduced. 
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radar  velocity  information,  but  the  disturbance  of  the  system  by  the  forward  loop  is  not 


compensated  for.  Therefore  Kx  must  be  made  zero,  if  the  system  shall  be  unsensitive 
against  external  disturbances.  But  by  omitting  Kx  the  whole  principle  of  the  simple 
platform  has  been  abandoned,  because  the  second  order  platform  with  84.4  minutes 
period  requires  components  of  the  highest  accuracy. 

4.  Combined  System  with  Short  Period  2nd  Order  Platform 

Since  one  could  not  obtain  a  correct  transient  response  with  an  essentially  first 
order  system  (i.e.  a  system  considerably  influenced  by  ),  the  question  arises,  if  one 
would  have  better  chances  with  a  non-tuned  short  period  second  order  platform.  Indeed, 
this  is  the  case.  Also  in  this  case  the  overall  system,  due  to  the  cross-feeds,  as  one 
loop  which  must  be  tuned  to  84-4  minutes  period.  As  is  well  known,  this  is  the  indis¬ 
pensable  requirement  for  each  system  with  correct  transient  response. 

The  system  is  shown  in  Figures  9  and  10.  The  additional  loops  are  inclosed  in  a 
dotted  box.  If  everything  contained  in  the  dotted  box  is  omitted,  Firgures  9  and  10  becomes 
identical  with  Figures  6  and  8.  The  complete  system  represented  in  Figure  10  has  the 
advantage  that  the  disturbing  influence  of  on  the  transient  response  can  be  compensated 
for  by  the  proper  choice  of  the  other  constants  (gains). 

The  transfer  function  can  be  deduced  in  the  same  way  as  it  was  done  in  paragraph  III. 2. 
Thus  one  finds  for  a  system  with  the  timing 


k2  =  w“ 

c2=  k3 

c3=  1- 


(35) 
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the  equation 


„  u  U2  +  KiK*s  +  Ki) Ci  (EDo)V  +  s  {s  +  KJC,}  (E)a 

V  =  V  +  — - ^ - - - - - - -  (36) 

s9  +  (Ci  +  KiK3)  s2  +  (CJCJC.  +  W2)  s  +  (CJCx  +  WjKJC  J 

The  transfer  function  for  the  angle  of  the  misalignment  could  also  be  deduced,  but  it 
shall  be  omitted  here;  however,  it  may  be  mentioned  that  the  main  term  becomes  equal 
to  zero,  if  equations  (35)  are  fulfilled.  That  is  evident,  because  equation  (36)  can 
be  valid  only  if  the  misalignment  of  the  platform  is  equal  to  zero  (except  for  the  errors). 

The  overall  system  described  by  (36)  will  be  stable  if 

Ki  <  +  KiKafCx  +  Ka)  (37) 

This  condition  imposes  certain  restrictions  on  the  selection  of  system  parameters.  If,  for 
instance,  a  short  period  inertial  platform  shall  be  used  then  equation  (37)  requires  that 
K  3  (a  minimum  platform  damping  term)  have  the  form 


This  is  required  to  render  overall  system  stability.  Thus,  positive  as  well  as  negative 
platform  damping  will  make  the  combined  system  stable.  However,  the  negative  values 
for  K  3  have  no  real  significance  because  it  would  be  very  impractical  to  use  a  highly 
suitable  inertial  component  in  the  overall  system. 

A  final  remark  shall  be  made  concerning  the  required  accuracy  of  the  system.  It 
is  evident  that,  in  order  to  make  the  84.4  minutes  loop  operate  correctly,  one  needs 
very  high  grade  components.  If  however  the  components. are  of  a  somewhat  lower  grade, 
the  accuracy  of  the  system  deteriorates  somewhat.  But  this  shows  up  only  in  the  error 
term  of  (36),  the  main  transfer  function  remains  essentially  equal  to  unity. 


(38) 
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jCHAPTER  IV 
(APPENDIX) 

MODE  OF  OPERATION  AND  ERRORS  OF  THE  DOPPLER  RADAR 

The  speed  measurement  with  the  doppler  radar  AN/APN-66  designed  by  General 
Precision  Laboratory  is  based  on  the  doppler  frequency  shift  of  radar  pulses  reflec¬ 
ted  from  the  ground.  The  device  has  an  array  of  four  antennae,  two  of  which  are  look¬ 
ing  right  and  left  forward,  and  die  other  two  are  looking  right  and  left  aft.  The  forward 
and  aft  frequencies  are  combined  crosswise.  In  this  way  one  obtains  and  error  signal 
if  the  antenna  array  is  not  aligned  with  the  direction  of  the  ground  path.  That  error 
signal  is  used  to  servo  the  antenna  array  into  the  direction  of  the  ground  path. 

Thus,  the  ground  velocity  vector  is  indicated  in  po(ar  coordinates,  i.e.  the  doppler 
radar  gives  the  magnitude  (ground  speed)  and  the  angle  of  the  ground  velocity  vector 
against  airplane  structure,  or  against  any  other  reference  direction  which  can  be  un¬ 
iquely  related  to  airplane  structure. 

The  speed  indication  of  the  doppler  radar  may  be  written: 

(Do)  V(B  _  A/C)  -  V(E  m  A/C)  +  (EDo)  .  A/C )  (t)  (39) 

where  A/C)  Speed  of  the  aircraft  against 

earth  (true  ground  speed) 

(EDo)  \f(E  m  a/C)  (t)  Error  of  doppler  speed 
indication 

The  error  is  com  posed  of  a  constant  bias  and  a  random  fluctuation  (noise).  The  bias  is 
claimed  to  be  very  small  if  a  careful  calibration  is  made  by  using  a  sufficiently  long  dis¬ 
tance  to  average  out  the  fluctuations.  However,  if  the  aircraft  is  flying  over  water,  the 
indication  of  the  ground  speed  will  be  falsified  by  ocean  currents  and  may  be  influenced 
by  wave  motions. 
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The  doppler  noise  is  caused  by  the  vibrations  of  the  antenna  array  and  by  the 
changes  in  terrain  and  the  consequent  changes  in  the  reflection  of  the  earth's 
surface.  If  (edo)  V  (t)  is  a  record  of  the  doppler  noise  taken  over  a  time  interval 
-T  -  t  -  +  T  then  the  Fourier  transform  of  this  record  can  be  written  as: 


(EDo)V(w)  =  3  [  (edo)V(t)  ]  = 


+  T 


if  (edo)  V  (t)  e  2  **  “\h 
-T 


(40) 


The  spectral  density  of  the  noise  is  defined  as: 

G,__  *>>  =  lim  “  I  (ED°)  V  (w)  : I2 
(EDo)V  I-®  J 


(41) 


Equation  (41)  can  be  approximated  sufficiently  well  by  taking  T  large  and  omitting  the 
limiting  process. 

Based  on  the  spectral  density  (41)  the  RMS  —  error  of  velocity  indication  can  be 
computed.  If  it  is  assumed  that  (  sdo)  V  (t)  is  not  cross-correlated  to  V(t),  the  in¬ 
fluence  of  the  doppler  noise  can  be  investigated  independently  from  the  dynamic  system 
response.  Designating  the  v, lighting  function  for  the  noise  (i.e.  the  filter  response 
function  for  unit  impulse  input)  by  W(t)  and  defining^ 


t 


<E>V, „,,(•> 


( edo)V(t)  •  W(t  -t)  di 


o 

the  noise  transfer  function  Y  (s)  is  found  by  taking  the  Laplace  transform  of  (42)  and 
applying  the  real  convolution  theorem: 


(42) 
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oc 


Y(s)  =  £  [W(t)  ]  =  W(t)e  Stdt  =  —  — 

J  (EDo)V(s) 


The  steady  state  component  of  the  noise  transfer  function  is  obtained  from  equation  (43) 
by  putting  S  =  2  Tti  (i),  he.  by  substituting  the  Fourier  transform  for  the  Laplace  transform. 
Therefore: 


Y (2  mi  w)  =  3  [W(t)  ]  =  W(t)e  2i;,Wtdt 


The  spectral  density  of  the  velocity  indication  error  is  then  found  to  be: 


°(E |Y(2tiu)l  'G(EDo)V(“> 

And  the  RMS  —  error  in  V,  ,  can  finally  be  computed  from: 


-  G 


(E)Vf 


(w)  dw 


|  Y  ( 2  tx i  to)  |  ’G(EDo)y(w)d(4  <*6> 


A  similar  derivation  holds  for  the  statistical  error  of  the  angle  (  edo)  i//(t)  of  the 
ground  velocity  vector  against  airplane  structure  (also  called  the  drift  or  crab  angle). 

Since  noise  records  for  (sdo)V(t)  and  (edo)  */>(t)  were  not  available  to  the  authors, 
no  numerical  computations  could  be  made.  However,  some  conclusions  regarding  the 
RMS  —  error  in  V can  be  drawn  from  theoretical  results  obtained  by  F.  B0  fierger 
in  M  Random  Errors  in  Doppler  Systems,  M  reference  4,  where  a  smoothing  time  of 
18  sec  was  found  to  be  sufficient  to  reduce  the  RMS  —  error  in  doppler  speed  to  less 
than  0.1%  of  true  ground  speed.  Also  the  RMS  —  error  in  true  ground  heading  was 


computed  to  be  less  than  0.1  degrees  after  6  sec  smoothing  time.  (Both  results  were 
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obtained  for  V  =  200  mph).  Since  the  inertial  platforms  usually  have  much  longer 
time  constants  or  periods,  resp.,  it  is  anticipated  that  they  will  act  as  very  ade¬ 
quate  filters  for  the  two  components  (edo)V(t)  and  (edo)  ^(t)  of  the  doppler 
indication,  so  that  the  resulting  RMS  —  error  in  is  expected  to  be  much 

less  than  0.1  %  in  magnitude  and  0.1  degree  in  direction.  For  the  8 4.4  minute 
system  the  doppler  noise  should  be  of  no  importance. 
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